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1. Introduction 

Any dimensionally -regularized |jl]] multiloop Feynman diagram with propagators 1 / {p^ — m^) 
can be written in the form of a finite sum of multiple Mellin-Barnes integrals ^ |3|] obtained via a 
Feynman-parameter or "a" representation: [Q] 



F {ajs,bkm,Ci,dj,x) = I dz\...dzr J ' , Vx -^l ^' ■ ■ --^r (1-1) 



where ajs^b^m G Q,c,-, S C, and 7 is chosen such that the integral exists. In this integral, the Xa 
are rational functions of kinematic invariants (masses and momenta) and the matrices ajs^b^r and 
vectors c,,<iy depend linearly on the dimension n of space-time.'. In dimensional regularization, n 
is treated as an arbitrary parameter, so that « = / — 2£ for a positive integer /. 

Formally,^ this integral can be expressed in terms of a sum of residues of the integrated ex- 
pression 



F [ap,bkm,c,d,a,xj =Y,Bax"^{f,d;x) , (1.2) 

a 

where the components of the vector a are defined in terms of the matrix components ajs^bkm and 
vectors c,d, the coefficients are ratios of F-functions with arguments depending on a, and the 
functions <I>(7; a;f) have the form 

withju«fa,v«/, gQ, 7,-,a;tGC. 

Let ?y = (0, • • • , 0, 1 , 0, • • • , 0) denote the unit vector with unity in its entry, and define = 

mi 
^1 



y,"' • • - a;'"' for any integer multi-index m = (mi , • • • , m^). In accordance with the definition of Refs. 



\ [\5i |16| ], the (formal) multiple series T.'^^qC {m)x'" , is called hypergeometric if, for each / = 1, • • • ,r. 



'starting from the a-representation for Feynman diagrams with (irreducible) numerators |^], the Melhn-Bames 
representation can be written exphcitly in terms of Symanzik polynomials.!^] In Refs. |^ it was shown that this 
a-representation can be understood as a linear combination of scalar Feynman diagrams of the original type with 
shifted powers of propagators and space-time dimension multiplying a tensor factor depending on external momenta. 
Through fifty years of the evaluation of Feynman diagrams, many auxiliary programs have been created for the genera- 
tion of Symanzik polynomials. In particular, the matrix representation for these polynomials derived by Nakanishi (see 
Eqs. (3.23)-(3.34) in Ref. |^) has been realized on FORM by Oleg Tarasov in Ref. [|], by Oleg Veretin in and 
recently in Ref . ||l^. 

^This is true under the condition that there is a common domain of convergence. For particular values of 
the kinematic variables and powers of the propagators, the Mellin-Bames representation may contain terms like 
T{a — s) X r(a-l-s) or r^(a — i) xT^{b s), where a and b are parameters and s is an integration variable. Two al- 
gorithms for the practical construction of the e-expansion of Mellin-Bames integrals with such singularities have been 
described in Refs. ||ll|], These are now implemented in several packages [|l2|. It will be quite interesting to apply modem 
mathematical algorithms, as was done in Ref. to the problem of singularities in the Mellin-Bames representation 
as well. In the framework of our approach, such singularities can be regularized by introducing an additional analytical 
regularization for each propagator (see Ref. [|l4|]) or by introducing masses to regulate IR singularities. General theo- 
rems on the properties of Feynman diagrams in dimensional regularization guarantee that a smooth limit of the auxiliary 
regularization should exist. 
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the ratio C(m + e,) /C{m) is a rational function in the index of summation (mi , • • • , m, ). Ore and 



Sato 017[ ] found (see also Ref. []16|]) that the coefficients of such a series have the general form 



where N,M > 0, Ay, 5y, 7y G C are arbitrary complex numbers, Hj, Vt : Z'' ^ Z are arbitrary integer- 
valued linear maps, and R is an arbitrary rational function. A series of this form is called a Horn- 
type hypergeometric series. 



We deduce from Eqs. (1.3) and (1.4) that, assuming there is a region of variables where each 
multiple series of Eq. (1.2) is convergent, any Feynman diagram can be understood as a special case 
of a Horn-type hypergeometric function with the functions R{m) equal to unity.^ One interesting 
property of Horn-type hypergeometric functions is the existence of a set of differential contiguous 
relations between functions with shifted arguments. We consider such relations in the following 
section. 

2. Contiguous Relations via Linear Differential Operators 

Let us consider a formal series of type ([0|). The sequences 7 = (/i , • " ' > 7/f ) ^rid a = (ai , • • • , Oj,) 
are called upper and lower parameters of the hypergeometric function, respectively. Two functions 



of type ( 1 . 3 ) with sets of parameters shifted by a unit, ^{Y+ec;d;x) and ( 7; a ; x) , are related by 



a linear differential operator: 



^{Y+e-;a-x)= y " 'y"j=^^ .^c=..j --. x'f--.x^ 



mi,m2,---,m,-— 
00 

E I 

mi,m2,--,m,-=0 



iLi + Yc) nf^.r (I^^i At;,m, + Yj) 



X^ ■ --x"' 



= l^caXa-^ + Ycj ^{f,o;x) = i7+^^^^jj<I>(7, a,x) . (2.1) 
Similar relations also exist for the lower parameters: 

'nJ^ir(i:tiAt7«'«« + 7; 



^{Y-d-ec;x) = LCvcfo/Wfo + a, 



mi,m2,...,m,=0fo=l \^k=l^ (111=1 Vktmb + Ok) 

Y VcbXh-^ + Oc\^(f,d;x) =L^^^^^^ ^^^{f,d;x) . (2.2) 



The linear differential operators U^^^ ^j, _j are called step-up and step-down operators for 

the upper and lower index, respectively. If additional step-down and step-up operators U^^, L+ 
satisfying 



Using the technique presented in Ref. ||l8(|, this statement can also be shown to be valid for the phase-space integral. 
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(i.e., the inverses of U^, L^ ), are consttucted, we can combine these operators to shift the param- 
eters of the hypergeometric function by any integer. This process of applying ?7^,L^^ to shift the 



parameters by an integer is called differential reduction of a hypergeometric function ( |1.3[ ). 

Algebraic relations between functions 0(7, a;x) with parameters shifted by integers are called 
contiguous relations. The development of systematic techniques for obtaining a complete set of 
contiguous relations has a long story. It was started by Gauss, who found the differential reduction 



for the 2F\ hypergeometric function in 1823. [ ]19| ] Numerous papers have since been published [ |20| ] 
on this problem. An algorithmic solution was found by Takayama in Ref. [p]], and those methods 



have been extended in a later series of publications [23|. 

Let us recall that any hypergeometric function can be considered to be the solution of a proper 
system of partial differential equations (PDEs). In particular, for a Horn-type hypergeometric func- 
tion, the system of PDEs can be derived from the coefficients of the series. In this case, the ratio of 
two coefficients can be presented as a ratio of two polynomials. 



C(fh + ej) _ Pj{m) 
C{m) ~ Qj{m) ' 

so that the Horn-type hypergeometric function satisfies the following system of equations: 



(2.3) 



= Dj{f,d,x)(i>{f,d,x) 



<I>(7,a,x), j = l,...,r. 

(2.4) 

Let 9^ be the left ideal of the ring T) of differential operators generated by the system of 
differential equations for a hypergeometric function (2A), Dj{f,a,x),j = l,.---,r. The first step 



in Takayama's algorithm is the construction of a Grobner basis (5 = {G/(7, d,x),i = 1, • • • of 
^R. Then the step-up operator corresponding to and step-down operator corresponding to L^^ 
are solutions to the linear equations 



TI+ - 1 



£QGK7,a,x) + [/| 

(=1 

££,G,(7,a,x)+L+^_j^^^]L[-^^_^^_i, = 1 



(2.5) 
(2.6) 



!=1 



where Q,^,- are arbitrary functions. This system has a solution if the left ideal generated by 6 U 
{?/+} is equal to D (see details in Ref. [^). 

In this way, the Horn-type structure provides an opportunity to reduce hypergeometric func- 
tions to a set of basis functions with parameters differing from the original values by integer shifts: 



Po{x)^{f+t,d + l;x) = £ P„,i,...m,( 



0(7; 6-; x) 



(2.7) 



dxi J \ dxr 

where Pq{x) and Pmi,-,m (x) are polynomials with respect to 7, a and x and k, I are lists of integers. 



The problem also can be solved via an Ore algebra [g^l approach to the relevant system of linear differential and 
difference (shift) operators. 
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3. Differential Reduction in Practice 



In real physical problems, the variables are generally not linearly independent: some of them 
can be equal to one another Xi = Xj,i ^ j. In this case, the Horn-type hypergeometric function 
generally is not expressible in terms of Horn hypergeometric functions with fewer variables. An- 
other important physical case is when variables belong to the surfaces where coefficient Po{x) of 
Eq. (p77|) vanishes (in the one- variable case it corresponds to z = 1). In all of these cases, the dif- 



ferential reduction cannot be directly applied. But if the l.h.s. of Eq. (2.7) is defined for that limit, 



the smooth limit of the r.h.s. of Eq. ( |2.7| ) will exist too. 

The problem is then to find this smooth limit. Here, "physics" plays a role. For the evaluation 
of physical processes, exact results in terms of hypergeometric function are not necessary, but only 
the coefficients of a Laurent expansion around an integer value of the space-time dimension. From 
that point of view, to guarantee that smooth limit exists, it is enough to prove that a smooth limit 
exists for all coefficients of the all-order e expansion. 

Recently, physicists have proven several theorems on the all-order e expansion of hypergeo- 



metric functions about integer and/or rational values of parameters ||24|, ^ |2^, gVp. A remarkable 
property of this construction is that for special values of parameters, the coefficients are expressible 
in terms of multiple polylogarithms [^]. For functions of this type, the limiting procedure is well 
understood. Unfortunately, existing theorems are not adequate to cover all values of parameters for 
hypergeometric functions generated by Feynman diagrams [27l 29]. 



4. Discussion and Conclusions 

We have presented formal arguments that any Feynman diagram can be treated as a finite sum 
of Horn-type hypergeometric functions. This applies, in particular, for off-shell diagrams with 
different values of masses in internal lines. In the physically interesting cases when some of the 
arguments are equal to one another, or belong to some singularity surface, a limiting procedure 
should be constructed. To find this limit, and strongly prove that the proper limit exists, the all- 
order £ expansion of hypergeometric function around rational values of parameters can be used. 

The Horn-type hypergeometric functions possess useful properties: the system of differential 
equation they satisfy is enough (i) for reduction of original function to a restricted set of basis 
functions (the number of basis functions follows directly from the system of equations); (ii) for 
the construction of the all-order e expansion for the basis hypergeometric functions in form of 
the Lappo-Danilevky solution. The first part of this algorithm has been discussed in [30, the 



second part in [26, 27]. The validity of our approach has been confirmed (in particular cases) by 



full agreement with the evaluation of the first coefficients of the e expansion constructed in [B3], 



by theorems about the all-order e expansion of hypergeometric functions proven in [25] with the 



help of another technique, and by comparisons of the results of the differential reduction of some 



Feynman diagrams with the results of reductions obtained via computer programs [33]. 

The above-mentioned properties of the hypergeometric representation, in particular Takayama's 
reduction algorithm, demonstrate that the hypergeometric representation is a universal tool for the 
evaluation of Feynman diagrams. Obtaining a decomposition of the integration region into regions 
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where each individual term is well-defined is one of the main problems in the hypergeometric 
approach (beyond the one-loop and one-fold cases). 

Series representations (in four dimensions) have been used in Ref. [34] to obtain a system 
of partial differential equations for A'^-point one-loop diagrams. In Ref. [^5|], the hypergeometric 
representation has been derived for an A/^-point one-loop diagram with arbitrary powers of the prop- 
agators. The possibility of using a hypergeometric representation for the reduction of Feynman 
diagrams was considered in Ref. [30]. The idea of using the Grobner basis technique for the re- 
duction of Feynman diagrams has been proposed by Tarasov [|3^ and received further extension in 



Ref. [37]. For practical applications of the "limiting" procedure, the first few coefficients of the e 
expansion are necessary. The e expansion is implemented in several packages [^] for a restricted 
class of hypergeometric functions. Some results for the finite harmonic sums are available in [39]. 

Another important class of developments includes techniques such as integration-by-parts 
[40], generalized recurrence relations and the differential equation approach [41], which make 
it possible to work directly with parameters of the Feynman diagram (the l.h.s. of Eq. (1^)) with- 
out splitting the diagram into a linear combination of Horn-type hypergeometric functions. Such 
direct analyses of Feynman diagrams as hypergeometric functions are based on the properties of 
dimensional regularization [|l|] , which treats all types of singularities (IR and UV) simultaneously. 

Some partial results of our research have been used in Ref. [|4^. A more detailed discussion 
will be presented in a forthcoming publication [43]. 
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